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Abstract 

This study is a numerical study on the coagulation as well as the sedimentation effect of nanofluids using the 
Brownian dynamics method. Three cases are simulated, focusing on the effects of the sizes, volume fraction, and 
£ potentials of nano-particles on the formation of coagulation and sedimentation of nanofluids. The rms fluctuation 
of the particle number concentration, as well as the flatness factor of it, is employed to study the formation and 
variation of the coagulation process. The results indicate a superposition of coagulation and sedimentation effect 
of small nano-particles. Moreover, it is stable of nanofluids with the volume fraction of particles below the limit of 
"resolution" of the fluids. In addition, the effect of £ potentials is against the formation of coagulation and positive 
to the stability of nanofluids. 



Introduction 

The nanofluid is characterized by the fluid with nan- 
ometer-sized solid particles dispersed in solution [1], 
which can increase the heat transfer coefficient [2-6], 
enhance the critical heat flux in boiling heat transfer 
[7-9], reduce the wall friction force [10], improve the 
optical characteristics [11], etc. Nano-sized particles are 
utilized because of its better stability than the suspen- 
sion of micro-sized particles. For a badly stable suspen- 
sion, sedimentation or coagulation (agglomeration) may 
occur. It compromises the above-mentioned advantages 
of the nano-suspension. 

As is well known [12], the occurrences of coagulation 
and sedimentation are the two main factors for the 
instability of nanofluid. The phenomenon of coagulation 
is characterized by the formation of particle clusters, i.e., 
particles are in contact with each other and the cohe- 
sion takes place. Then, the clusters grow up. Many 
researchers investigated the coagulation effect of parti- 
cles by the Brownian dynamics simulation, focusing on 
the formation of gelation [13], coagulation rates [14], 
particle network [15], etc. For example, Hutter [14] 
identified the characteristic coagulation time scales in 
colloidal suspensions, and measured their dependencies 
on the solid content and potential interaction para- 
meters. He also deduced different cluster-cluster 
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bonding mechanisms in the presence of an energy bar- 
rier, etc. Besides, the sedimentation always occurs after 
a big particle cluster is established, i.e., the particles 
within the cluster sediment flow downward because of 
the increased effect of the gravity of the cluster over the 
buoyancy force of it, and reduced the effect of Brownian 
motion to the big cluster. Many researches were devoted 
to the sedimentation [16-18] using the Brownian 
dynamics simulation too. For example, Soppe and 
Jannsen [17] studied the sediment formation of colloidal 
particle by a process of irreversible single-particle accre- 
tion. They used the algorithm of Ermak and McCammon, 
incorporating the inter-particle forces and hydrodynamic 
interaction on the two-particle level, and analyzed the 
effect of two-particle hydrodynamic interactions on the 
sediment structure, etc. They found that the process of 
sediment formation by colloidal particle is the result of a 
delicate balance of sediment field strength, DLVO inter- 
actions, and hydrodynamic interactions. 

However, there is an important issue about which few 
researches have been concerned: the interaction between 
the coagulation and sedimentation for the instability of 
nanofluids. For example, the processing of coagulation 
causes the particle clusters to grow up, and then the 
large clusters are more prone to sedimentation than that 
of small clusters because of the intensive gravity effect. 
In other words, the coagulation effect is able to augment 
the sedimentation effect. Thus, this study is intended to 
carry out some research on this issue, exploring the 
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complex interaction as well as the close relation 
between the coagulation and the sedimentation 
phenomena. 

Numerical model 

Governing equation 

In this study, the Brownian dynamics technique is 
employed to investigate the motion of nanoparticles. 
The governing equation is the so-called Langevin equa- 
tion, which is formulated as follows [19]: 
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where the superscript 0 indicates that the variable is 
corresponding to the beginning of the time step At; r t is 
the /th component of the position vector of particle; 
is the element of the diffusion tensor indexed by (i, j); Fj 
is the force experienced by the yth particle; /c B is the 
Boltzmann constant; and T denotes the temperature. 
The displacement Ri(At) is a random displacement of a 
Gaussian distribution with a zero expectation and a 

( 2D -At ) covariance ( (R^A^R^At)) = 2D? At ). In this 

study, the Rotne-Prager tensor [20] is utilized as approx- 
imations to the hydrodynamic interaction: 
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where r\ is the viscosity, a is the particle radius, 3% is the 
Kronecker delta, is the vector from the center of parti- 
cle i to the center of particle and J is the unit tensor. 

Moreover, three forces, i.e., the attractive Van der 
Walls force, f v , the repulsive electrostatic force by the 
electric double layer, f e , and the gravity force f g are con- 
sidered as the forces experienced by any particle, which 
are formulated as follows [12]: 
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where A, d, e, k, £ p f , p, and g are the Hamaker con- 
stant, the particle diameter, the electric permittivity of 
the fluid, the inverse of the double-layer thickness, the 
zeta potential of the suspension, the density of the fluid, 
the density of the particle, and the gravity acceleration, 
respectively. 

It is noted in Equation (4) that the results for r t j - d = 

0 is meaningless when the contact between the two par- 
ticles occur, and they will adhere to each other or 
rebound back. Thus, we treat the condition with 

1 = ^-1 — — < T as the situation when the two particles 

d 

are separated, so that Equation (4) works. Otherwise, it 
results in coalescence between the two particles. Once 
the coalescence between colliding particles takes place, 
the clusters start growing up. 

Simulation conditions 

In this study, three cases with different diameters of par- 
ticles (Case 1), different volume fractions (Case 2), and 
different zeta potentials (Case 3) are simulated, respec- 
tively (Table 1). For these cases, the parameters of the 
material, as well as other relevant parameters, are illu- 
strated in Table 2. 

In this simulation, the boundary conditions in the x 
and y directions (Table 2) in the horizontal plane are 
both periodic, whereas the top and bottom walls of the 
simulation domain in the z-direction are treated as adhe- 
sive walls to which the particles adhere immediately once 
they come into contact with them. It is reasonable to 
conclude thus, since the agglomerated particle clusters 
always adhere to the bottom walls or the top interfaces. 

Table 1 Three cases with different diameters of particles, 
volume fractions, and zeta potentials 

Case 1: under different diameters d 



Case 



Case 



a) d 0 

b) d, 

c) d 2 -- 



10 nm, N p 
■ 25 nm, N p 
50 nm, N n 



1200, i//= 0.153, C = 0.0 eV 
: 1200, y/= 0.153, £ = 0.0 eV 
1200, (//= 0.153, ^ = 0.0 eV 



2: under different volume fractions yj 

a) d y = 25 nm, N p = 400, y/ = 0.051, £ = 0.0 eV 

b) d } = 25 nm, N p = 1200, y/ = 0.153, 5 = 0.0 eV 

c) d } = 25 nm, N p = 2100, y/ = 0.268, £ = 0.0 eV 

d) d ] = 25 nm, N p = 4200, yi = 0.537, ? = 0.0 eV 
3: under different zeta potentials £ 

a) d ] = 25 nm, N p = 1200, yj = 0.153, £ = 0 eV 

b) d } = 25 nm, N p = 1200, y/ = 0.153, C, = 0.01 eV 

c) d } = 25 nm, N p = 1200, y/ = 0.153, £ = 0.025 eV 

d) d } = 25 nm, N p = 1200, y/ = 0.153, £ = 0.05 eV 
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Table 2 Parameters used in this simulation 



Simulation domains {L Xl L y , L z ) 
Mesh sizes {3 X , S y , 3 Z ) 
Temperature T (°C) 
Hamaker constant A (J) 
Diameters of particle {d 0l d h d 2 ), (nm) 
Density of fluid p f (kg/m 3 ) 
Density of particle (p p ) (kg/m 3 ) 
Viscosity of fluid r\ (Pa s) 
Simulation time step At (ns) 



(10, 25, 50) 
0.993 x 10 3 
6.4 x 10 3 



25 



2 x 10" 



{Ad, Ad, 256d) 
{Ad, Ad, Ad) 



1.0 x 10" 



10 



r3 



Initially, a random distribution is given to the parti- 
cles. As time advances, the possible movements of parti- 
cles are computed through solution of the governing 
equations. 

Simulation results 

Case 1: effect of particle sizes 

This section deals with the effect of particle sizes on the 
coagulation and sedimentation. Figure 1 shows the simula- 
tion results at t = 0, 5, 10, and 50 [is for Case 1 due to the 
effects of different sizes of particles. Figure la, b, c, d shows 
the results of d = 10 nm at t = 0, 5, 10, and 50 (is, respec- 
tively. Similarly, Figure le, f, g, h shows the results for d = 
25 nm, whereas Figure li, j, k, 1 shows them for d = 50 nm. 
It is seen that the coagulation takes place the most 
intensively and rapidly for the smallest size of particles (Fig- 
ure la, b, c, d), moderately for the intermediate size of par- 
ticles (Figure le, f, g, h), and weakly and slowly for the 
largest size of particles (Figure le, f, g, h, i, j, k, 1). More 
importantly, the results of the intermediate sizes are due to 
coagulation but with weak sedimentation, whereas the 
results of the smallest sizes are due to both the effects of 
coagulation and sedimentation. For the largest size of the 
particles, it is neither due to coagulation nor sedimentation. 
It looks complicated. As is known, the larger particles bear 
the major effect of gravity, and they are the most prone to 
sedimentation. However, it is only true of the single particle 
without coagulation. With the superposition of the coagula- 
tion effect, it can amplify or augment the trend of sedimen- 
tation through coagulation. Owing to the increasing 
agglomeration of the particles, the gravity effect may play 
an important and even a dominant role, which causes pos- 
sible sedimentation of the whole agglomeration (the upper 
part of the agglomeration in Figure Id is due to the adhe- 
sive boundary on the upper wall). In other words, there 
exists a balance between the sedimentation effect of the 
large-sized individual particles and the sedimentation effect 
of small-sized aggregated particles. The former is caused 
solely by the gravity effect, whereas the latter is caused by 
the superposition of the coagulation and the gravity effects, 
i.e., the amplification and augmentation of the gravity 
effects of the aggregated particles due to coagulation. 
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It is necessary to mention that Figure li, j, k, 1 does 
not indicate the stability of the nanofluids. Alternatively, 
it indicates a relatively stable characteristic compared to 
Figure la, b, c, d. After the evolution over a long time, 
possible coagulation or sedimentation may also occur. 

In order quantify the degree of coagulation, we need 
to define some functions. Let us divide the simulation 
domain L x x L y x L z by N x x N y x N z cubic meshes by 
the cell volume (S x x S y x S z ). The mean number con- 
centration c of particles is the mean number of parti- 
cles within each mesh volume (S x x S y x d z ). Then, the 
rms value of the particle concentration Ri and the flat- 
ness factor of the number concentration R 4 are formu- 
lated as follows: 



R x = f-^- J"(c(x,y,z) - c) 2 dz 



V/2 



1 



-j- J(c(x, y, z) - c) 4 dz 



(7) 



(8) 



The rms of concentration means the fluctuation of the 
number concentration of particles, and it is closely related 
to the formation of particle clusters due to coagulation. 
The flatness factor means the intensity of fluctuation of 
the number concentration, thereby indicating the intensity 
of coagulation. Thus, these two functions are helpful in 
enabling the quantification of particle coagulation. 

Figure 2a, b shows the Ri and R 4 for Case 1. It is seen 
from Figure 2a, b that the coagulation of the small parti- 
cles is the fastest. They become almost totally coagulated 
immediately even at the beginning. Comparatively, the 
coagulation of the larger particles takes place slowly and 
increasing steadily. However, the final level of coagula- 
tion of the larger particles is greater than that of smaller 
particles. In addition, the degree as well as the rapidity of 
the coagulation of the intermediate particles is intermedi- 
ate between that of the smaller and the larger particles. 

Case 2: effects of volume fractions 

In this section, the effect of volume fraction, i.e., the 
concentration of particles, is studied. As aforemen- 
tioned, the smaller particles are more prone to coagulate 
than the larger particles, under the same condition of 
the volume fractions. However, the process of coagula- 
tion is also closely related to the number of particles 
contained in it. 

For example, Figure 3a, b, c, d shows one of the 
results of Case 2 where only n p - 400 particles are 
simulated. Compared to Figure la, b, c, d, it is seen that 
the coagulation does not takes place at t = 50 [is. It says 
that the coagulation is to be regarded appropriately as a 
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Figure 2 The flatness factor of the concentration distribution 
of nanoparticles. (The R^a) and R 4 (b) for Case 1). 



process with an excess of particle content. When the 
particle numbers go beyond the superior limit of the 
resolvent, then the coagulation will certainly take place. 
Thus, when the n p = 4200 particles are simulated, more 
intensive coagulations are observed correspondingly 
(Figure 3e, f, g, h). 

In addition, Figure 4a, b shows the Ri and R 4 of Case 2. 
It is seen from Figure 4 that the Ri and R 4 for n p - 400 
are always relatively of small value, indicating a stable 
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Figure 3 Snapshots of simulation results (Case 2 for n p = 400 
(a-d) and n p = 4200 (e-h) at t = 0, 5, 10, and 50 \as respectively). 



status almost without coagulation and sedimentation, 
although R 4 is slightly fluctuated when t < 0.06. More- 
over, compared to n p = 1200, it is seen that the concen- 
tration fluctuation R x and flatness factor of concentration 
R 4 for n p = 400 are relatively of lower values. It validates 
the conclusions derived from the observation of Figure 
3a, b, c, d, in comparison with Figure le, f, g, h. 

With the increased number of particles, it is seen that 
the Ri and R 4 are increased too (n p - 1200 and 2100, 
respectively, Figure 4). However, when the particle num- 
ber is extremely large, all the spaces are almost stuffed 
with particles, leading to a homogeneous distribution 
and a low fluctuation in the number concentration 
(rip = 4200, in Figure 4). 
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Figure 4 The flatness factor of the concentration distribution 
of nanoparticles. (The R^a) and R 4 (b) for Case 2). 



Case 3: effects of £ potentials 

The previous sections showed the results with £ = 0 eV. 
As seen from Equation (5), no repulsive effect has been 
considered between the particles since f e = 0. Thus, this 
section will focus on the effect of the repulsive effect by 
varying the £ potentials. 

Comparing with Figure le, f, g, h, it is seen that the 
degree of coagulation with £ = 0.01 eV is attenuated 
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(Figure 5a, b, c, d), and it almost disappears with £ = 0.05 
eV (Figure 5e, f, g, h). It indicates that the repulsive effect 
induced by the £ potentials is beneficial to the stability of 
nanofluids, since it acts against the coagulation process. 

This conclusion is also validated by the variations of 
Ri and R 4 (Figure 6a, b, respectively). It is seen that the 
fluctuation of the number concentration of particles 
with £ = 0.05 eV increases much more slowly than the 
cases with smaller Although the R x and R 4 for £ = 
0.05 eV are still increasing with time, which indicates 
that the coagulation process may still take place after a 
fairly long time, their effects against the formation of 
coagulation are very clear. In other words, it is positively 
beneficial for the stability of nanofluids. 



3 



Figure 5 Snapshots of simulation results. (Case 3 for £ = 0.01 eV 
(a-d) and £ = 0.05 eV (e-h) at t = 0, 5, 10, and 50 pis, respectively). 
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Figure 6 The flatness factor of the concentration distribution 
of nanoparticles. (The R^a) and R 4 (b) for Case 3). 



Conclusions 

The findings of this study are briefly summed up as follows: 

1. A complicated superposition of the coagulation 
and sedimentation effects for small particle is 
observed. The mechanisms of sedimentation for the 
larger and the smaller particles are different. The 
former is caused mainly by the great gravity effect of 
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any individual particle, whereas the latter is mainly 
due to the coagulation process, and the superposi- 
tion of coagulation causes the sedimentation of the 
whole agglomeration of particles. 

2. There exists a superior limit of the fluid for parti- 
cle content. When the volume fraction is below the 
limit, it is hard for the coagulation to occur. In con- 
trast, the coagulation will certainly take place when 
the concentration of nanoparticles is beyond the 
capacity of "resolution" of the fluids. 

3. The effect of £ potentials is beneficial for the sta- 
bility of nanofluid, since it resists the formation of 
coagulation. In other words, increase in the value of 
£ potentials is helpful to make the nanofluid more 
stable. 
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